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Abstract 

We study experimentally and theoretically the influence of light polarization on the interference pat- 
terns behind a diffracting grating. Different states of polarization and configurations are been con- 
sidered. The experiments are analyzed in terms of electromagnetic energy (EME) flow lines, which 
can be eventually identified with the paths followed by photons. This gives rise to a novel trajectory 
interpretation of the Arago-Fresnel laws for polarized light, which we compare with interpretations 
based on the concept of "which- way" (or "which-slit" ) information. 

Keywords: Arago-Fresnel laws, which-way information, electromagnetic-energy-flow line, photon inter- 
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1 Introduction 



Experiments with polarized light are well known since the beginning of the 19th century. By this time, 
Arago and Fresnel [1] found a series of results related to Young's interference experiment when the latter 
is carried out with polarized light. These results were summarized in the form of four laws, known as the 
Arago-Fresnel laws [1-4]. According to these laws, two beams of the same linear polarization interfere with 
each other just as natural rays do. However, no interference pattern is observable if the two interfering 
beams are linearly polarized in orthogonal directions. In the 1970s, a series of experiments with linearly 
and elliptically polarized laser light were performed, illustrating the Arago-Fresnel laws [5,6]. Actually, 
the original Arago-Fresnel laws were generalized to include elliptically polarized light: two beams with 
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the same polarization state interfere with each other just as natural rays do, but no interference pattern 
will be observable if the two interfering beams are elliptically polarized, with opposite handedness and 
mutually orthogonal major axes. 

The standard interpretation given to the disappearance of interference after inserting mutually or- 
thogonal polarizers after the slits is usually based on the Copenhagen notion of the external observer's 
knowledge (information) about the photon paths, i.e., the slit traversed by the photon in its way to 
the detection screen. If, after allowing the two diffracted beams to reach the screen, additional polar- 
izers are inserted between the slits and the screen, interference emerges again [5,7]. This reappearance 
is attributed to the erasure of the observer's which-way information, in this case actually "which-slit" 
information [8,9]. 

Recently, lecture demonstrations of interference [10] and quantum erasure [11,12] with single photons 
have been experimentally realized by Dimitrova and Weis. In both cases, experiments were first carried 
out with strong light. Then, in a second step, the same experiments were performed with strongly 
attenuated light by means of single photon counting. In this way the formation of the corresponding 
pattern by the progressive detection of single photons was observed, as it has also be done with matter 
particles [13,14]. Motivated by the feasibility of this kind of experiments, Sanz et al. [15] have recently 
proposed a description based on electromagnetic energy (EME) flow lines to analyze the formation or 
disappearance of interference features depending on polarization. As shown, within this theoretical 
framework it is possible to understand experiments like the one just mentioned on the basis of the 
topology displayed by the EME flow lines in their way from the grating to the screen and on how the 
presence and arrangement of polarizers behind the grating influences them. This approach, suggested as 
a tool to analyze and explore basic experiments in quantum optics, thus combines the ideas and methods 
from the hydrodynamical formulation of Maxwell's equations [16] and Schrodinger's wave mechanics. A 
similar synthesis has also been exploited by M. Man'ko and coworkers to study charged particle beam 
optics [17] and to describe the classical propagation of electromagnetic waves through optical fibers [18]. 

Earlier, EME flow lines for a particular polarization of light were determined by Braunbek and 
Laukien [19, 20] in the case of a half plate or, later, following a similar approach, by Prosser [21-23] 
for the case of a single and double-slit interference. More recently, some of us [24] considered the case 
of linearly polarized light and provided a method to systematically compute the evolution of the EME 
flow lines by means of the transverse momentum approach [25]. Also, Gondran and Gondran [26] have 
considered the EME-flow-line approach to study diffraction by a circular aperture and a circular opaque 
disk, the latter giving rise to the well-known Poisson-Arago spot phenomenon. Authors put this result in 
the historical perspective by showing that EME flow lines provide a complementary answer to Frenel's 
answer to the questions presented in 1818 by the French Academy "deduce by mathematical induction 
the movements of the rays during their crossing near the bodies" [27]. 

The four Arago-Fresnel laws governing the interference of polarized light were determined experimen- 
tally using natural and linearly polarized light. As mentioned above, experimental veriflcations of these 
laws were done with intense laser light. Here, we analyze these laws under the conditions of single-photon 
count experiments, providing a description and explanation in terms of EME flow lines. Accordingly, 
the organization of this work which encompasses both experiment and theory is as follows. Thus, the 
experimental setup and findings are described in Sec. [2j In Sec. [3l in order to be self-contained, we 
present a brief overview of the EME-flow-line approach, described in detail in [15]. In particular, we 
focus on how the presence of polarizes influences both electromagnetic (EM) fleld and EME flow lines. 
In Sect. |4]we show the EME flow lines obtained with the conditions used in the experiments. Finally, in 
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Figure 1: Schemes of the experimental setups used. Abbreviations: P: linear polarizer, DS: double-slit, 
BC: Babinet compensator, DPS: diffusive projection screen, RSBP: rotating slit beam profiler. 



Sec. [5] we provide a discussion on the interpretation of the experimental results in terms of EME flow lines 
as well as a comparison between this interpretation and that one based on the concepts of "which-slit" 
information and quantum erasure [8,9]. 



2 Influence of the polarization state on the double-slit interference 
pattern: Experimental facts 

Experiments were performed at the University of Fribourg. Light from a green laser module (A = 
532.5 nm) was sent through a polarizing system consisting of a linear polarizer and a Babinet compensator 
(BC), which allowed the realization of arbitrary linear and elliptical polarizations (Figs. [1^ and b). The 
beam polarization was measured with a repositionable commercial polarimeter that was removed for the 
actual diffraction measurements. The double-slit consisted of two 18 mm long slits of 100 ^m, whose 
centers were separated by 250 ^m. For the experiments shown in Fig. [21 the laser beam has a nearly 
circular Gaussian intensity profile, 

I{r) = Ioe-^^'/^\ (1) 

with w = lA mm. 

Two-dimensional interference patterns were recorded with the setup of Fig. [1^ for various incident 
light polarizations. The diffraction patterns were projected onto a thin sheet of paper and their intensity 
distribution photographed by a digital camera. The results are shown in Fig. [2] for eight circular, elliptical 
and linear polarization states. In a second experiment, the interference pattern produced by the expanded 
beam was recorded by a scanning slit beam profiler (Fig. [H)) for eight polarization states. The signal of 
the beam profiler corresponds to the x-dependence of the intensity integrated over the z-direction. 

The results of Fig. [2] seem to indicate a slight polarization dependent transversal shift. This could 
be traced back to a systematic beam displacement during adjustments of the Babinet compensator. The 
results of Fig. [3]were obtained after resolving this problem. Fig. [2] — and more convincingly Fig. [3] — show 
that the interference pattern does not depend on the state of polarization of the incident laser light. In 
a third experiment — carried out with incident circularly polarized light — , we have inserted a horizontal 
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Figure 3: Interference patterns from a double-slit (0.1 mm slit width, separated by 0.25 mm) along the 
x-direction [in mm] obtained from the intensity (EME density) integrated over the z-direction for eight 
polarization states (shown as insets) of the incident light. 
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Figure 4: Intensity (EME density) as a function of x at the screen behind the double-slit (0.1 mm slit 
width, separated by 0.25 mm) illuminated with circularly polarized light in two cases: with H and V 
polarizers (blue curve) just behind each slit and without them (purple curve). 



and a vertical polarizer after each slit, respectively. Figure H] shows the diffraction patterns obtained with 
and without these polarizers. The peak intensities are in the ratio 4:1, as expected. We note that single 
photon interference experiments with mutually orthogonally polarized waves are discussed in [11,12] in 
the frame of quantum erasure processes. 
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3 Theoretical background 



3.1 EM field and EME density beiiind a double-slit grating for different incident 
polarizations 

Consider the EM field behind a grating situated on the X Z-plane, at y = 0, with shts parallel to the 
z-axis, whose width (5) is much larger along the z-direction than along the x-direction. Accordingly, 
we can assume the EME density to be independent of the z-coordinate and, therefore, the electric and 
magnetic fields will not depend either on this coordinate. As shown in the literature [20], this allows 
us to express the electric and magnetic fields as a sum of S-polarized (-Be.x = -E'e,?/ = -f^e,^ = 0) and 
i?-polarized {Hfi,x = Hh,y = Eh^z = 0) components, i.e., 

E = Ee + H,, = Ee + — [V X H,,] , 

ueo (2) 

H = He + H;, = [V X Ee] + H,,. 

Let us assume that the incident EM field is a plane wave of frequency u and wave vector k propagating 
along the y-axis, with arbitrary polarization. It can be written as [15] 



Eo(r,t) = [-/3e*('=S'+'^)x + ae^'^f z] e 
Ho(r,t) = 



-iu)t 
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(3) 



where a, (3 and (j) are real quantities. When (/> = or vr and a and /3 are arbitrary, the polarization of 
the incident wave is linear. If = ib7r/2 and a = /3, the polarization is circular. In all other cases the 
polarization is elliptic. In ^ we recognize the ii^-polarized and i?-polarized components, 

Eo,e = a^'oz, 
V A'o 

with 

^o(r) = e*'=^ (5) 

being the corresponding scalar field satisfying Helmholtz's equation. 

If the boundary conditions at the grating for £'e^2(r) and Hh^zi^) are the same and both fields satisfy 
Helmholtz's equation, we may also assume that these fields are proportional to a scalar field, ^'(r), which 
also satisfies the same boundary conditions and the Helmholtz equation. Accordingly, we can express 
these fields as 

Ee = a^z, 

(6) 



V /^o 



Within the paraxial approximation and for completely transparent slits with their support being totally 
absorbing, the solution to Helmholtz's equation behind a double-slit reads as 



^{x,y) = Mx,y) = \h^ e-^-'^e^^y [ ,0-)e''''^'-^'^"/''dx' , (7) 
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where ^'0(2;', 0~) is the incident scalar wave just before the grating. Taking into account 
we find that 

E = El + E2, H = Hi + H2, 

where 



H, 



k dy k dx 

ia dtpi ^ ia dihi ^ 

y + TT- y + 



y + aViZ, 
kf3e''^ 



and dZD, 
(8) 

(9) 



ufiQ dy uj^Q dx oj^o 
represent, respectively, the electric and magnetic fields propagating from the ith slit, with = 1,2. 

As seen above, experimentally one measures the light intensity, i.e., the EME density on the screen 
positioned at some distance from the diffraction grating. From the general expression for the EME energy 
density, 



U{r) 



1 



[eoE(r).E*(r)+^oH(r).H*(r)], 



(10) 



we find that the EME density of the incident EM wave ([3]) is constant (i.e., independent of both x and y), 



(11) 



However, as a result of the interaction between the EM field and the grating, the EME density behind 
this grating becomes dependent on x and y, 



The solution ^ to the Helmholtz equation satisfies the approximate relations 
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(12) 



dx 



d^ 

dy 



d"^ 
dy 



ik^ 



and, therefore, the EME density (|12p becomes proportional to i.e.. 



C/(r) = |(a2 + /32)|^|2^|(„2 + ^2) 



|^'iP + |^'2p + 2Re(V'iV2) 



(13) 



(14) 



As can be noticed from this expression, the x and y dependence of the EME density can be determined 
through \^\^- The polarization of the incident EM wave has no influence on the dependence of the 
EME density on x and y and, therefore, the interference pattern will be independent of the incident 
field polarization. This distribution is the same for linear, circular or elliptic polarization, as seen in the 
experimental results displayed in Figs. [2] and [3j It is important to stress that this conclusion arises from 
assuming that the EM field does not depend on the z-coordinate. 

The independence of the EME density on polarization implies that the interference pattern will be the 
same for both linearly polarized and natural light, in agreement with one of the empirical findings observed 
by Arago and Fresnel, namely the first Arago-Fresnel law. This law states that "two rays polarized in 
one and the same plane act on or interfere with each other just as natural rays, so that the phenomena 
of interference in the two species of light are absolutely the same" [1,3]. The generalization of this law 
also includes elliptic polarization, and circular polarization as a special case of elliptic polarization [4]. 
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3.2 EM field and EME density behind a double-sHt grating followed by orthogonal 
polarizers 

The approach described above can also be apphed to the study of the EM field and EME density behind 
two slits which are followed by two linear orthogonal polarizers. Now, instead of considering Eqs. ([6]), 
the E'-polarized and ff-polarized components of the EM field are expressed [15] in terms of ipi and ip2, 
as 




^0 iS I - 



(15) 



Thus, by substituting ()15p into ([2]), we obtain the EM field behind the two slits covered by the orthogonal 
polarizers, 



Ej = — -— X — y + aipiT., 

k oy k ox 

Hj = — y H — y H ip2'^- 

w/xq oy w/xq ox W/Xq 

From this EM field, the expression for the EME density will read as 



(16) 



Uir) 



Co 2 

— a 
2 



2 



a 



dx 



+ 



dy 



d'ip2 
dx 



+ 



dip2 
dy 



(17) 



Under the presence of orthogonal polarizers behind the slits, the EME density becomes the simple sum of 
EME densities, which spread out independently from each slit. As a consequence, interference fringes are 
absent. This is in agreement with the second law experimentally established by Arago and Fresnel [1], 
namely the second Arago- Fresnel law [3]. This law states that "two rays primitively polarized in opposite 
planes (i.e., at right angles to each other) have no appreciable action on each other, in the very same 
circumstances where rays of natural light would interfere so as to destroy each other" . The generalization 
of this law also includes elliptic and circular polarization of opposite handedness [4]. 



4 EME-flow-line description 

The properties of the EME density behind a grating and the generalized Arago-Fresnel laws may be 
better understood with the aid of EME flow lines, which can be interpreted as photon trajectories [15]. 
The equation from which the EME flow lines are obtained is 

dr _ 1 S(r) 

d^~^U{^' ^^^^ 
where S(r) is the real part of the complex Poynting vector, 

S(r) = iRe[E(r)xH*(r)], (19) 
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U{r) is the time-averaged EME density and s is the certain arc-length of the corresponding path, 
substitution of Eqs. dZ])-© into ([T9]) renders the trajectory equations along each direction, 



The 



dx 
ds 

dy 
ds 
dz 
ds 



4:kU{r) 

AkU{r) 

ieoa/3 sin cj) 
2A;2C/(r) 



dx 

^^^)* 

dy 
0^ dy^f* 



dx 
dy 

d^* 



(20) 



dx dy dy dx 



The projection of the EME flow lines on the Xy-plane behind a double-slit grating illuminated by 
circularly polarized laser light, determined from Eqs. (j20p and the corresponding experimental conditions, 
are plotted in Fig. [5l As can be noticed, the EME flow lines (photon trajectories) accumulate according 
to U{r) [15], as given by Eq. ()14p . It is important to stress that, in the case of circular and elliptic 
polarization, EME flow lines do not remain within the XY-p\a,ne, but there is a flow of energy parcels 
(photons) along the vertical axis. Nevertheless, the projection of these three-dimensional trajectories onto 
the Xy-plane does not depend on the phase angle, (/>. Due to this fact, the resulting distribution of end 
points associated with the trajectories along the x-axis does not depend on polarization. This is consistent 
with the conclusions derived in the previous section from analysis of the EME density dependence on 
polarization. 

The number of bright and dark fringes obtained in the experiment, where d = 0.25 mm, 5 = 0.1 mm 
and A = 532.5 nm, agrees fairly well with the corresponding number predicted by the probability ampli- 
tude of transverse momenta. 



c{kx 



1 



dx'^'(x',0+)e-^'=-^' 



1 



2^ ^=l,2■'A 



dx'^o{x' ,0-)e 



(21) 



which, in the case of two slits, is given by 



c(/cx) 



2 sin(/c^5/2) 



Tr6 



cos{kxS/2). 



(22) 



The centers for the zeroth (central), first (±1) and second (±2) bright fringes are determined through 
the relation cos{kxS/2) = ±1, which renders kxfi = 0, /cx,±i = ±27r/d and kx^±2 = ±47r/d, respectively. 
On the other hand, the centers of the dark fringes are obtained taking into account that cos{kx6/2) = 
and, therefore, kx = =b(2n + l)Tr/d, with n = 0, 1, 2, . . .. 

The intensity of the bright fringes decreases because of the envelope of the interference pattern, which 
is given by first factor and describes the single-slit diffraction. This factor vanishes as kx increases from 
zero to the value determined by the relation sin(A;^5/2), i.e., kx = ib27r/5. The point where the intensity 
falls to zero due to diffraction will coincide with one of the interference minima if 



2tt _ {2n + l)7r 
T ~ d 



d _2n + l 
'5 ~ d ' 



(23) 



As seen in the experimental results, this condition holds, because d/5 = 5/2 = 2.5. Therefore, the second 
bright fringe will be very weak (almost fainting) and the third dark fringe will be negligible in comparison 
with the remaining continuous dark pattern, as seen in Figs. [2] and [31 
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Figure 5: Top: EME density at L = 558 mm from a double-slit grating illuminated by a circularly 
polarized laser beam of wavelength A = 532.5 nm. Centre: 30 EME flow lines behind a double-slit 
grating determined from (I20p and the scalar solution ([7]) of Helmholtz's equation. Bottom: Enlargement 
of the EME flow lines in the near field region. The slits are assumed to be completely transparent, with 
their support being completely absorbing. The distance between the center of the slits is d = 0.25 mm 
and the slit width is 6 = 0.1 mm. 
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Figure 6: Top: EME density at L = 558 mm from a double-slit grating with orthogonal polarizers upon 
the slits and illuminated by a circularly polarized laser beam of wavelength A = 532.5 nm. Centre: 
30 EME flow lines behind a double-slit grating, determined from Eqs. (I16p . (llSp and (jlOh . and the 
scalar solution ([7|) of Helmholtz's equation. Bottom: Enlargement of the EME flow lines in the near 
field region. The slits are assumed to be completely transparent, with their support being completely 
absorbing. The distance between the center of the slits is d = 0.25 mm and the slit width is (5 = 0.1 mm. 
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When mutually orthogonal polarizers are positioned upon the slits, both the EME distribution and 
the EME flow lines change dramatically [15]. The resultant distribution of trajectory end points along 
the X-axis will be consistent with expression (jl7p for the EME density and interference fringes will be 
absent. 

5 Summary and conclusions: EME flow lines versus "which-slit" in- 
formation 

The analysis presented above throws some novel light on the first and second Arago-Fresnel laws. More- 
over, it also results relevant to the discussion on the wave-particle duality of light. The fact that fringes 
disappear after covering the slits with orthogonal polarizers has been used by the proponents of the 
principle of complementarity to affirm that information about the path destroys the interference [8,9]. 
Within this type of argumentations or interpretations, "path" just means the slit through which the 
photon passes and any discussion about the whole path, from the grating to the screen, is thus forbidden. 

Argumentations based on "which-slit" information are closely related to the reasoning leading to the 
notion of quantum erasers [8,9]. By applying this reasoning to the sequence in which the first experiment is 
a double-slit experiment with orthogonal polarizers and the second experiment is a double-slit experiment 
without orthogonal polarizers, one could assign the lack of interference to the existence of "which-slit" 
information. The appearance of interference after removing the polarizers could be then attributed to 
erasing and, consequently, to the loss of information about the slit which the energy parcel or photon 
passed through. 

The interpretation used here considers a set of entire paths from the grating to the screen, such paths 
being determined from the EM field and the Poynting vector. EME flow lines starting from slit 1 will end 
up in the left-hand side part of the screen, while those starting from slit 2 will end up in the right-hand 
side, as also shown in quantum mechanics for matter particles [28]. This is valid both when interference 
is present and when no interference fringes are observed [29,30]. But, the distribution of this EME 
flow lines is different in the two cases. In the absence of polarizers, the distribution shows interference 
fringes (Fig. [5]); in the presence of polarizers, the fringes are absent (Fig. [6|), in perfect agreement with 
experimental results shown at Fig. [H Hence, observer's information about the slit through which photon 
(energy parcel) went through is not relevant for the existence of interference. What is relevant it is the 
form of the EME fleld and its corresponding form, which will model consequently the distribution of 
trajectories and their topology. 
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